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摘 要 : 本 文 以 单 级 蜗 壳 式 离 心 泵 为 研究 对 象 ， 基 于 RANS 进行 了 全 流 场 的 非 定常 数值 计算 , 探讨 了 叶轮 的 非 定 常 
流动 特征 ， 重 点 分 析 了 叶轮 出 口 沿 轴 向 变化 的 圆周 面 流域 内 各 节点 的 压力 分 布 特性 。 结 果 表明 : 在 叶片 压力 面 及 叶 
片 后 缘 尾 迹 区 发 生 了 剧烈 的 满 流 脉动 现象 ， 不同 工 况 下 叶轮 出 口 圆 周 方向 的 压力 系数 G, 均 在 隔 舌 附近 产生 极 值 ; 
额定 工 况 Q/QN=1L 时 ， 叶 轮 出 口 圆周 方向 的 压力 系数 C, 呈 稳 定 周期 性 分 布 ， 随 着 工 况 的 改变 ， 压 力 系 数 C, 沿 蜗 壳 
周 向 分 布 不 均匀 性 增加 ; 小 流量 工 况 下 叶轮 流域 内 发 生 了 明显 的 流动 分 离 现象 ， 流 道内 产生 了 较 大 尺度 的 分 离 涡 
破坏 了 叶轮 出 口 圆周 方向 压力 系数 C, 的 周期 性 分 布 特性 ， 叶 轮 出 口 圆周 压力 系数 C, 沿 中 截面 两 侧 的 流域 内 几乎 对 
称 分 布 ; 在 Q/QN=0.2 小 流量 工 况 下 ， 压 力 脉动 频谱 图 低频 段 内 出 现 了 较 多 复杂 的 激励 信号 ， 认 为 这 与 流域 内 的 分 
离 涡 结 构 具 有 一 定 关 联 。 
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Study on the Internal Flow Characteristics of the Centrifugal 
Pump under Non-design Conditions 
GAO Bo WANG Xing-Ning LIZhong ZHANG Ning YANG Min-Guan 
(School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China) 

Abstract: Set the single stage volute centrifugal pump as research object meanwhile based on the RANS to 
simulate the unsteady flow characteristics of the whole flow field, and then researched on the unsteady flow 
characteristics of the impeller, which focused on the pressure distribution and fluctuation characteristics of 
each node along the axial direction that changes in circumference field of impeller outlet. The results shows 
that: the phenomenon of intense turbulence pulsation happened in the pressure surface and the trailing edge of 
the blade; the pressure coefficient C, in the circumferential direction of the impeller outlet reached to the 
extreme value around the tongue under different conditions; the pressure coefficient C in the circumferential 
direction of the impeller outlet has a stable periodic distribution under the rated condition, the non-uniformity 
of the pressure coefficient C, along the circumferential direction increases while the working condition 
changed; the flow separation phenomenon occurred obviously under small flow rate condition and appeared 
large scale separation vortex which destroyed the periodic distribution of the impeller outlet pressure 
coefficient C, along the circumferential direction. Complicated exciting signals were found in the pressure 
pulsation frequency graph at small flow rate condition, which may related to the vortex structure in the pump. 
Key words: Centrifugal pump; unsteady flow; unsteady pressure distribution; flow separation 


站 | 阁 影响 并 导致 机 械 部 件 疫 劳 过 度 ， 引 起 振动 和 噪声 。 

因此 ， 深 入 分 析 离 心 泵 内 部 的 压力 分 布 及 脉动 特性 
离心 泵 作为 一 种 基本 的 能 量 转换 和 流体 运输 效 将 成 为 迫切 关注 的 焦点 59。 

置 广泛 应 用 于 工农 业 ， 航 空 推进 和 日 常生 活 中 吓 。 不 国内 外 学 者 对 离心 泵 内 部 不 稳定 流动 进行 了 | 

对 称 结构 的 蜗 壳 和 高 速 旋转 的 叶轮 ， 在 流体 间 产 生 泛 研 究 ， 表 建 平 等 中] 采用 LES 数值 模拟 了 离心 泵 内 

了 动静 干涉 作用 ， 导致 离心 和 的 运行 个 稳定 记性 及 六 部 的 非 定常 流动 结构 ， 认 为 内 部 压力 脉动 变化 明显 

生 随 时 间 变化 的 压力 脉动 5 ， 这 将 对 泵 的 性 能 产生 是 由 于 动静 干涉 所 致 ，Wangl 等 基于 LES 模型 对 离 
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心 和 叶片 进口 及 蜗 壳 基 圆 线 的 压力 脉动 特性 进行 了 
数值 研究 及 试验 验证 ， 沿 蜗 壳 基 圆 各 监测 点 的 压力 
分 布 周 期 与 叶片 数 一 致 。 总 的 来 说 ， 由 于 离心 泵 内 
部 的 非 稳 态 流 场 的 复杂 性 ， 其 内 在 机 理 仍 处 于 计算 
探索 和 试验 研究 阶段 。 

本 文 基于 RANS 并 采用 一 @ 满 流 模 型 对 单 级 蜗 
壳 式 离心 泵 的 内 部 流动 进行 了 三 维 非 定常 的 数值 计 
算 及 分 析 ， 重 点 探讨 了 叶轮 出 口 的 非 定常 流动 特性 ， 
同时 采用 统计 法 分 析 了 叶轮 出 口 沿 轴 向 变化 其 圆周 
面 内 各 节点 的 压力 分 布 特性 。 通 过 不 同 工 况 下 的 对 
比分 析 ， 揭 示 了 离心 条 叶轮 出 口 流域 内 的 非 定常 流 
动 特性 ， 其 结果 可 为 离心 汞 内 部 非 定常 流动 的 研究 
提供 一 定 的 理论 支撑 。 


1 参数 及 数值 计算 方法 
1.1 “模型 泵 参数 


文中 所 用 离心 泵 基本 性 能 参数 如 表 1 所 示 ， 叶 
轮 和 蜗 壳 的 三 维 实体 均 采用 Pro/E5.0 造型 。 
表 1. 离 心 泵 基本 参数 


Tablel Basic parameters of centrifugal pump 


主要 参数 数值 
设计 流量 QN (mh) 29 
扬程 H/m 9.8 
转速 n (r/min) 1450 
比 转速 ns 86 
叶片 数 Z 6 

叶轮 出 口 宽度 bymm 14 
叶轮 进口 直径 Dj/mm 80 
叶轮 出 口 直径 Dymm 184 


1.2 ”数值 计算 方法 
1.2.1 网 格 划分 及 计算 模型 


整个 模型 由 一 个 动 叶轮 水 体 、 一 个 静止 蜗 壳 水 
体 及 进出 口 延 伸 段 组 成 ， 并 将 上 述 四 个 子 域 整体 作 
为 数值 计算 求解 区 域 。 

本 文采 用 ICEM 网 格 划分 软件 对 整个 计算 域 进 
行 了 结构 化 网 格 划分 ， 整 个 流 道 区 域 所 划分 网 格 总 
数 为 1128281, 其 中 叶轮 部 分 网 格 数 为 566716, 蜗 壳 
部 分 网 格 为 435128。 经 网 格 无 关 性 检验 可 知 ， 当 网 
格 数 大 于 110 万 以 后 , 随 着 网 格 数 的 增多 , 数值 模拟 
的 结果 显示 扬程 和 效率 波动 性 很 小 。 因 此 ， 计 算 采 
用 的 物理 模型 及 离散 格式 具备 了 空间 收敛 性 和 数值 
稳定 性 ， 达 到 了 计算 的 精度 要 求 。 数 值 计 算 区 域 的 
局 部 网 格 划分 如 图 1 所 示 ， 为 了 得 到 更 为 精准 的 数 
值 计 算 结 果 ， 分 别 在 流动 区 域 径 向 截面 所 处 位 置 如 
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叶片 周围 及 蜗 壳 隔 舌 等 位 置 进行 了 网 格 的 质 
密 ， 一 定 程度 上 调整 了 网 格 的 节点 数 。 


了 加 


大 Sa pa 
图 1 模型 泵 网 格 示意 图 
Figure.1 Mesh of the model pump 


1.2.2 ”灌流 模型 的 选择 


为 了 研究 不 同 注 流 模型 对 计算 结果 的 影响 ， 本 
文 分 别 采 用 标准 ke 模型、k 一 & 模 型 及 SST 模 型 对 
相同 边界 条 件 下 泵 的 能 量 性 能 进行 计算 ， 并 与 模型 
泵 的 实验 结果 进行 了 比较 分 析 。 图 2 为 模型 泵 在 不 
同 工 况 下 采用 四 种 湛 流 模型 计算 所 得 结果 和 实验 结 
果 进 行 对 比 ， 其 中 w=gH/u2? 表示 无 量 纲 扬 程 系 数 ， 
横 坐 标 为 运行 工 况 点 流量 与 额定 流量 的 比值 ， 由 
Q/QN 表 示 。 


图 2 模型 泵 性 能 曲线 


Figure.2 Performance characteristics of pump 
通过 对 比 发 现 ， 所 采用 的 三 种 满 流 模型 计算 所 
得 能 量 性 能 曲线 与 试验 结果 在 整体 变化 趋势 上 大 体 
致 。 当 运行 工 况 处 于 Q/QN=0.6 一 1.2 范围 内 时 ， 数 
值 计 算 扬程 值 和 试验 结果 吻合 较 好 ， 最 大 误差 小 于 
5%; 在 大 流量 工 况 下 ， 随 着 流量 的 增 大 ， 数 值 计算 
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误差 增 大 ， 但 通过 各 模型 所 得 计算 结果 之 间 却 相差 
不 大 ; 此 外 ， 由 效率 特性 曲线 可 以 看 出 ， 与 其 它 几 
ee 小 流量 工 况 下 k - w 模 型 所 得 计算 
结果 与 试验 结果 吻合 度 最 高 ， 考 虑 到 需 针 对 离心 泵 
外 特性 得 到 相对 准确 预测 吉 果 方面 的 需求 ， 并 综合 
计算 的 收敛 性 及 误差 等 因素 ， 最 终 汕 流 模 型 采用 标 
准 k - w 模 型 。 


2.2.3 ”边界 条 件 


通过 ANSYS CFX13.0 对 离心 泵 计算 域 进行 非 
定 芝 流动 计算 ， 采 用 多 重 坐标 系 求解。 壁面 采用 无 
滑 移 边界 条 件 ， 近 壁 区 域 采 用 标准 壁面 函数 处 理 
考虑 到 加 工 精 您 方面 的 有 限 性， 对 各 壁面 进行 粗糙 
度 定 义 为 0.Imm。 i 
口 设置 为 自由 出 流 。 在 进行 非 定 常 计算 时 ， 叶 轮 区 
域 相对 计算 域 其 它 网 格 运动 ， 将 动静 交界 面 设置 为 
瞬 态 动静 转子 模式 ， 这 对 于 计算 动静 干涉 作用 是 至 
关 重 要 的 。 非 定常 数值 模拟 将 定常 计算 结果 作为 初 
始 条 件 ， 取 叶轮 每 旋转 1° 作 为 一 个 时 间 步 长 ， 即 
AT=1.149x10“s， 每 计算 360 步 为 叶轮 旋转 的 1 个 周 
期 。 


3 计算 结果 分 析 


3.1 满 动 能 分 布 特性 
为 了 较为 直观 的 分 析 离 心 泵 内 部 产生 的 周期 性 


涡流 脉动 强度 ， 根 据 文献 中 定义 了 漠 流 强度 系数 Tu， 
计算 式 如 下 所 示 ;: 
天 9 和 
L010 = 0) 
- 
式 中 天 表示 滑动 能 ， 上 表示 叶轮 出 口 速度 。 


图 3 所 示 为 Q/Qw=0.2、1.4 工 况 下 离心 泵 叶轮 中 
间 截 面 处 某 时 刻 的 满 动 能 强度 系数 ZT 分布 图 ， 通 过 
坐标 对 叶轮 流域 展开 后 得 到 。 由 图 可 以 看 出 ， 在 
QQN=0.2 工 况 下 在 叶轮 流域 内 叶片 压力 面 附 近 的 汝 
动能 产生 较为 剧烈 的 变化 ， 且 叶片 后 缘 尾 迹 区 内 江 
动能 变化 较 大 ， 汕 流 不 稳定 性 明显 ， 此 时 发 生 了 剧 


cp 
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烈 的 满 流 脉动 现象 ， 能 量 耗 散 较 大 。 而 与 Q/QN=0.2 
工 况 相 比 较 后 可 知 ，Q/Qw=1.4 时 由 于 没有 发 生 流动 
分 离 因而 满 动 能 变化 不 大 。 

Tu hy 00 “ma 
0 260 se 
22.500 
18.750 
15.000 - 一 


| 11.250 

| 7.500 
3:750 
0.000 
图 3. Q/Qw=0.2、1.4 工 况 叶轮 流 道中 间 截 面 的 满 流 强度 系数 


Figure.3 Parameter distributions on middle plane of the impeller 


under Q/QN=0.2 and 1.4 


3.2 叶轮 出 口 的 压力 分 布 特性 


本 文 所 用 压力 分 布 采取 了 标准 化 形式 ， 即 定义 
了 无 量 纲 压 力 系数 。 通 过 计算 基于 叶轮 圆周 速度 UU， 
与 压力 标准 化 后 的 比值 得 到 所 求 无 量 纲 的 压力 系数 
C0,, 其 中 变量 Popz 轨 用 于 表示 计算 流域 内 任意 网 格 
节点 的 压力 ， 由 此 定义 的 无 量 纲 压力 系数 1 如 式 
(2) 所 示 : 


05pU? . 

为 了 得 到 叶轮 旋转 过 程 中 叶轮 出 口 圆 周 方向 节 
点 压力 随时 间 的 变化 规律 ， 通 过 全 流域 的 非 定 常 计 
算 ， 得 到 如 图 4 所 示 的 离心 泵 中 间 截 面 处 叶轮 出 口 
沿 轴 向 变化 的 圆周 面 流域 内 各 节点 的 压力 分 布 特性 
沿 圆周 方向 的 压力 系数 C, 分 布 。 由 图 可 以 看 出 在 离 
心 泵 运行 过 程 中 ， 不 同 工 况 下 的 叶轮 出 口 圆周 方向 
所 得 压力 系数 C, 均 在 隔 舌 附近 达到 了 极 小 值 ( 图 中 
出 现 极 小 值 点 为 9 = 26? 时 )， 通 过 分 析 可 知 这 与 隔 
天 处 强烈 的 动静 干涉 作用 相关 ;， 与 小 流量 工 况 下 所 
得 压力 系数 分 布 曲 线 相 比 , 在 额定 工 况 Q/QNw=1 时 叶 
轮 出 口 圆 周 方向 的 压力 分 布 呈 稳 定 的 周期 分 布 ， 其 
周期 与 叶片 数 一 致 ， 且 各 流 道 出 口 压力 分 布 均匀 ; 
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图 4. 不 同 工 
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圆周 方向 的 压力 系数 分 布 


Figure.4 Pressure parameter distributions on the cylindrical direction around the impeller under different working conditions 
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而 小 流量 工 况 下 ，0 > 26? 时 的 压力 系数 C, 沿 
圆周 方向 出 现 增 大 的 趋势 ， 且 Q/QN=0.2 工 况 下 更 有 
压力 分 布 周期 性 发 生 明 显 变化 ， 叶 轮 出 口 一 个 流 道 
内 产生 了 两 个 峰值 ， 这 是 因为 小 流量 工 况 下 流动 较 
为 复杂 ， 叶 轮流 域内 发 生 了 明显 的 流动 分 离 现象 ， 
流 道内 产生 了 涡流 所 致 ， 从 而 破坏 了 流动 的 稳定 性 
与 周期 性 。 

图 $ 所 示 为 与 图 4 工 况 相 对 应 的 叶轮 流域 内 的 
相对 流速 分 布 图 ， 可 以 看 出 在 小 流量 工 况 下 叶轮 流 
道内 出 现 了 大 尺度 的 流动 分 离 涡 ， 分 别 位 于 叶片 进 
口 、 叶 片 出 口 工作 面 及 叶片 出 口 背面 ， 部 分 流 道 流 
涡 明 显 严重 后 时 出 现 堵塞 流 道 现 象 ， 而 额定 工 况 时 
明显 看 到 相对 流速 分 布 均匀 。 
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本 文采 用 统计 法 分 析 了 叶轮 出 口 圆周 沿 叶 片 厚 
度 变化 的 流域 内 各 节点 的 压力 分 布 特性 ， 得 到 了 
QQN=1.4 工 况 时 压力 系数 治 叶轮 出 口 圆周 面 的 分 布 
如 图 6 所 示 。 从 图 中 可 以 看 出 ， 叶 轮 出 口 圆周 面 的 
压力 系数 Cr 在 隔 舌 前 后 具有 较 大 的 压力 梯度 ， 与 小 
流量 工 况 相反 ， 大 流量 工 况 叶 轮 出 口 圆周 方向 压力 
系数 Cr 在 隔 舌 前 表现 为 极 小 值 ， 在 隔 天 附近 出 现 了 
拐点 ， 紧 跟着 出 现 的 是 压力 系数 Co 的 极 大 值 ， 而 后 
压力 系数 C, 沿 圆周 方向 呈 递 减 趋势 ， 总 的 来 看 ， 压 
力 系 数 Cr 的 变化 周期 基本 可 划分 为 60"， 与 叶片 流 


道 数 一 致 。 叶 轮 出 口 圆周 压力 系数 C, 沿 中 截面 两 侧 
的 流域 内 几乎 对 称 分 布 。 


(a) Q/QN=0.2 


(b) Q/QN=0.6 
图 5. 不 同 工 况 下 叶轮 流域 内 相对 流速 分 布 


(c) Q/QN=1.0 


Figure.5 The relative velocity distribution in the impeller flow filed under different working conditions 
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图 6. Q/Qw=1.4 工 况 下 叶轮 


b 口 圆周 面 的 压力 系数 分 布 


Figure.6 Pressure parameter distributions on the cylindrical surface around the impeller at Q/QN=1.4 


3.3 蜗 壳 测 点 的 压力 脉动 频谱 分 析 


为 了 分 析 蜗 壳 内 部 由 于 内 流 场 的 不 稳定 流动 所 
引发 的 压力 脉动 特性 随 流量 的 变化 情况 ， 进 行 非 定 
常 流 场 计算 时 布置 监测 点 进行 时 域 信号 的 监测 ， 并 
最 终 通过 快速 傅 里 叶 变 换 得 到 对 应 的 压力 脉动 频 域 
图 如 图 7 所 示 ， 选 取 P4 和 P5 测 点 分 析 非 设计 工 况 
下 的 压力 脉动 特性 ,其 中 P4 测 点 表示 蜗 壳 第 IV 断 面 ， 
P5 测 点 表示 的 是 蜗 壳 第 工 断 面 。 频 域 图 中 横 轴 表示 
频率 值 f, 纵 轴 表 示 各 频率 下 对 应 的 压力 脉动 幅 值 A。 


本 文 所 用 叶轮 的 旋转 速度 n=1450r/min， 计 算 模型 泵 
的 轴 频 太 =24.2Hz, 且 叶 轮 叶 片 数 三 6, 故 fppp=145Hz。 
从 图 7 可 以 看 出 脉动 幅 值 明显 集中 在 低频 区 ， 压 力 
脉动 的 主峰 值 对 应 的 是 叶片 通过 频率 145Hz 处 ， 叶 
频 1 倍 频率 对 应 幅 值 明显 高 于 高 次 谐 波 下 幅 值 ， 即 
| 频 占 主导 作用 ， 且 同时 可 以 在 叶 频 整数 倍 处 发 现 
| 频 高 次 谐 波 处 均 对 应 着 一 定 峰 值 。 此 外 由 图 7 还 
可 看 出 ， 在 Q/QN=0.2 极 小 流量 下 ， 测 点 的 压力 脉动 
幅 值 较 Q/QN=1.4 工 况 较 大 ， 且 低频 段 内 出 现 了 更 多 


中 
中 


的 激励 频率 信号 ， 这 与 小 流量 工 况 下 内 流 场 的 复杂 
流动 特性 相关 ， 认 为 与 流域 内 的 分 离 涡 结构 上 共有 一 
定 关 联 。 
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(bb) QQN=1.4 工 况 下 蜗 壳 测 点 压力 脉动 频谱 
7. 不 同 工 况 下 P4、P5 测 点 压力 频谱 图 


Figure.7 Frequency domain Signals of different conditions 


Da 


distributions at point4 and point5 
4 结 论 


1) 不 同 工 况 下 叶轮 出 口 圆周 方向 的 压力 系数 Cy 均 在 
隔 舌 附近 产生 极 值 ， 在 额定 工 况 时 ， 叶 轮 出 口 圆周 
方向 的 压力 系数 C, 分 布 呈 稳定 周期 性 分 布 ， 且 各 流 
道 出 口 压力 分 布 均匀 ， 随 着 工 况 的 改变 ， 压 力 系数 
C; 沿 蜗 壳 周 向 分 布 不 均匀 性 增加 。 

2) 小 流量 工 况 下 叶轮 流域 内 发 生 了 明显 的 流动 分 离 
现象 ， 流 道内 产生 了 较 大 尺度 的 分 离 涡 ， 破 坏 了 叶 
轮 出 口 圆周 方向 压力 系数 C, 的 周期 性 分 布 特性 ; 

3 ) 叶 轮 出 口 圆周 压力 系数 C 在 隔 舌 前 后 位 置 具有 较 
大 的 压力 梯度 ， 且 C 沿 中 截面 两 侧 流 域内 几乎 对 称 
分 布 ; 
4) 在 Q/Qw-0.2 流量 下 ， 低 频段 内 出 现 了 更 多 复杂 
的 激励 频率 信号 ， 认 为 与 流域 内 的 分 离 涡 结构 具有 
一 定 关联 。 


也 | 
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